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Resumen 


La programación dinámica (PD) fue desarrollada por Richard E. Bellman en 
la década de los cincuenta (Bellman, 1957). Es una herramienta matemática que 
ha sido utilizada en campos tan diversos como ingeniería, inteligencia artificial, 
economía, etcétera (Dormido et al., 2002). El campo de la planificación de recursos 
hidráulicos no ha sido la excepción y se reportan aplicaciones con buenos resultados 
en diversos países. Tiene como limitante un alto costo computacional; necesita de 
grandes recursos tanto de tiempo como de memoria y la demanda de ellos crece 
exponencialmente, al incrementar el número de variables de estado. En este trabajo 
se presenta la forma en que se organizó un programa de cálculo para contender 
con este problema y obtener las políticas de operación Óptima de un sistema de tres 
presas que Operan en cascada, variando condiciones de penalización por déficit, 
derrame o por exceder las curvas guía (que son los niveles de almacenamiento que 
se desea no sean rebasados durante la operación para la seguridad del sistema y 
que establece la Comisión Nacional del Agua). Las políticas obtenidas se simularon 
usando tanto el registro histórico como registros sintéticos para probar su bondad 
en el funcionamiento del sistema. 


Palabras clave: programación dinámica, operación en cascada, políticas de 
Operación, curvas guía. 


Introducción 


La PD es una metodología matemática orien- 
tada a la solución de problemas en los que se 
deben tomar decisiones en un horizonte de 
planeación de N etapas sucesivas, teniendo 
como objetivo maximizar o minimizar el costo 
de esas decisiones (Ramos, 2010). Ha probado 
ser un método eficiente para determinar 
políticas óptimas de operación de sistemas de 
presas hidroeléctricas, en las que se concilia el 
propósito de maximizar la energía generada con 
el de minimizar las posibilidades de derrames 
o de déficit en el sistema. Así, en México se ha 
aplicado a los sistemas de presas del río Grijal- 
va (Domínguez y Mendoza, 2000; Domínguez 
et al., 2001a; Arganis, 2004; Domínguez et al., 


2006; Arganis y Domínguez, 2008; Domínguez y 
Arganis, 2009; Arganis et al., 2009), del río Fuerte 
(Domínguez et al., 2007), etcétera, y en países 
como Chile (Cámac et al., 1995), España (Ramos, 
2010; Estalrich et al., 2000), Perú (CESEL, S.A., 
2008), Colombia (Bedoya y Herrera, 2005), Brasil 
(Macedo y Franco, 2002), entre otros, se reportan 
aplicaciones exitosas de la PD a diversos siste- 
mas hidroeléctricos. Sin embargo, se considera 
que la aplicación del método está limitada por 
lo que ha dado en llamarse “la maldición del 
número de estados”, debido a que al crecer este 
número también lo hace el número de variables 
de decisión, el de barras del histograma de 
probabilidades de ingreso, etcétera, de tal 
forma que el número de cálculos se incrementa 
en una proporción exponencial. Para salvar 
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esta limitante se han desarrollado diversas 
variantes de la PD, tales como la dinámica 
incremental (ID) (Larson, 1968), la estocástica 
(Bertsekas, 1976) y la discreta diferencial (DD) 
(Jacobson y Maine, 1970; Heidari et al., 1971); 
cabe destacar que la DD es una generalización 
de la ID (Nopmongcol y Askew, 1976), y ambas 
requieren que la función objetivo cumpla con 
ser diferenciable y las restricciones impuestas 
deben ser lineales. 

En este trabajo se analiza un sistema de tres 
presas que funcionan en cascada; se requiere 
determinar políticas de operación que hagan 
óptimo su funcionamiento y se eligió la técnica 
de la programación dinámica en su versión 
estocástica para ello. El algoritmo se organizó 
de tal forma que se concilie la limitante del uso 
y demanda de los recursos computacionales. 

Una restricción adicional impuesta a las 
políticas que aquí se determinaron fue con- 
siderar que se debe reducir a un mínimo la 
probabilidad de que el almacenamiento en las 


presas del sistema supere el correspondiente a 
una curva denominada “curva guía”. 

Para comprobar la bondad de las políticas 
obtenidas, además de la simulación de su fun- 
cionamiento con el registro histórico y con 
registros sintéticos, se compararon los resul- 
tados que se obtienen con dos posibles 
volúmenes de discretización, a fin de verificar 
si la solución ya no mejora al reducirlo. 


Objetivo 


Desarrollar un algoritmo que permita obtener 
y comparar el comportamiento de distintas 
políticas de operación quincenal obtenidas 
utilizando la programación dinámica en su 
versión estocástica para el sistema de embalses 
del río Santiago, compuesto por tres presas 
que operan en cascada: La Yesca, El Cajón y 
Aguamilpa, considerando tanto las necesidades 
asociadas a la generación de energía como las 
correspondientes al control de inundaciones. 
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Figura 1. Sistema de presas del río Santiago (Jalisco y Nayarit). 
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Figura 2. Perfil de las presas hidroeléctricas del río Santiago (Jalisco y Nayarit). 


Elaborar escenarios considerando diversas 
opciones de curva guía para los niveles en 
dichas presas. 


Metodología 
Ubicación del sitio 


El río Santiago se localiza en el occidente de la 
república mexicana. Nace en el lago de Chapala, 
en el estado de Jalisco, a una altitud de 1 500 
msnm, y sigue su trayecto por los estados de 
Jalisco y Nayarit. A lo largo de este río se han 
construido presas, cuyo principal propósito es 
la generación de energía eléctrica. El sistema 
está formado por las presas Santa Rosa, La 
Yesca (actualmente en construcción), El Cajón 
y Aguamilpa. En la figura 1 se ilustra una vista 
en planta de la forma en que están dispuestas 
estas presas y en la figura 2 un perfil. 

En los cuadros 1 a 3 se presentan algunas 
de las características físicas más relevantes de 
las presas. 


Programación dinámica estocástica 
La determinación de las políticas de operación 


de un sistema de presas que opera en cascada 
es complicada; está ligada con procesos de 


Cuadro 1. Características físicas de la presa La Yesca. 


Dato Unidad Cantidad 
Volumen aproximado mó 12 300 000 
Embalse 
Elevación del NAMINO msnm 518 
Elevación del NAMO msnm 575 
Elevación del NAME msnm 578 
Capacidad útil x 10% m* 113592 
Elevación de la corona msnm 579 
Obras de generación 
Turbinas 
Número y tipo de unidades Francis 2 
Unidad de generación MW 375 
Gasto de diseño por unidad m*/s 153.6 
Generación media anual GWh 1210 
Elevación de desfogue m 390 


decisión secuencial que hacen necesario el uso 
de algoritmos de optimización para resolver 
el problema. El planteamiento de la función 
objetivo es trascendental, ya que debe conciliar, 
por una parte, los máximos beneficios por 
generación, y por la otra, establecer las pena- 
lizaciones por posibles condiciones de derrame 
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Cuadro 2. Características físicas de la presa El Cajón. 


Dato Unidad Cantidad 
Embalse 
NAME msnm 394 
NAMO msnm 391 
NAMINO msnm 346 
Capacidad útil x 10% m* 1 335.32 
Capacidad de control de avenidas x 10% m* 110.97 
Obra de generación 
Obra de toma Tipo rampa 
Número de tomas 2 
Nivel medio de desfogue msnm 221.24 
Turbinas 
Número de turbinas 2 
Tipo Francis, eje vertical 


Potencia nominal de la turbina a carga neta de 


os MW 375.00 MW 
Carga neta de diseño m 156.54 m 
Gasto de diseño por unidad m?/s 259.70 m*/s 
Cuadro 3. Características físicas de la presa Aguamilpa. 
Dato Unidad Cantidad 
Embalse 
NAME msnm 232 
NAMO msnm 220 
NAMINO msnm 190 
Capacidad útil x 10% m* 2 629 
Capacidad de control de avenidas x 10% m* 1 418 
Capacidad muerta x 10% m* 2 827 


Obra de generación 


Obra de toma 


SS Número de tomas 5 
a Factor de planta % 19.46 
uv 

9 Desfogue 
5 

5 Nivel medio de desfogue msnm 66.29 
S Nivel máximo de desfogue 

o Ñ pala msnm 69 
ys (tres unidades a carga máxima) 

5 Turbinas 

3 Número de turbinas 3 
=l Tipo Francis, vertical 
z Potencia nominal de la turbina MW 320 
El Carga neta de diseño m 145.1 
$ Gasto de diseño por unidad m/s 240 


C 


« vj30]0u72], 


Arganis-Juárez et al., Operación de tres presas hidroeléctricas usando curvas guía y programación dinámica estocástica 


o de déficit en el sistema. De manera adicional 
se debe considerar la naturaleza aleatoria 
de los volúmenes de ingreso a los embalses, 
incluyendo el grado de autocorrelación y de 
correlaciones cruzadas existentes en ellos. 

En este estudio se dividió el año en N 
etapas; se consideraron los tres embalses di- 
vididos en cierto número de estados, NS, 
dependiendo de la capacidad útil de cada uno 
y proponiendo un incremento de volumen 
AV para hacer la discretización del problema. 
Las variables de decisión o control son las 
extracciones, k, las cuales pueden tomar, para 
cada etapa, desde un valor mínimo hasta un 
valor máximo, dependiendo del volumen 
máximo turbinable en cada embalse; las 
variables aleatorias son los escurrimientos 
por cuenca propia, x, a las que se definió 
por medio de sus funciones de densidad 
de probabilidad, para obtener el beneficio 
esperado, que forma a la función objetivo del 
problema, las restricciones en el modelo están 
dadas por la ecuación de continuidad: j = 1 + 
x - k; con los estados 1 < 1¡<NS, las entradas 1 < 
x < nx; las extracciones: kmín < k < kmáx. 

La función objetivo que se propuso consis- 
te en la maximización del valor esperado 
del beneficio por generación, imponiendo 
penalizaciones en el caso de déficit, de de- 
rrame O bien de rebase de la curva guía 
propuesta para cada presa, quedando la 
siguiente expresión: 


NP 
FO = Máx EO) (EG ¡; -C jeyy Dert;; — C ge Def; =Cog,) 


1 
li=1 


(1) 


Siendo FO la función objetivo; NP, el número 
de presas (tres en este caso); ii, el subíndice 
para indicar la presa (í¡ =1 es La Yesca; ¡¡=2, El 
Cajón, ¡¡=3, Aguamilpa); E(), el valor esperado; 
EG, la energía generada en cada presa; C,,., 
Cay Y Cy los valores para los coeficientes de 
penalización por derrame, déficit y rebase de 


la curva guía en cada presa, respectivamente; 


Derr, la magnitud del derrame, y Def, la del 
déficit. 


Desarrollo del algoritmo y forma propuesta 
para resolverlo 


La política Óptima se obtiene usando las 
siguientes ecuaciones recursivas del método de 
la PD en su versión estocástica: 


a NS, NS, NS, 

Br (1, 1,13) = > x S A(bpar)+B,.1 (is 52,3) 
J=15,=1j3=1 

siendo: 


1=La Yesca; 2=El Cajón; 3=Aguamilpa 


A= nx rd 9,1, (o 72 M9 y,x, (13113) 
bpar = ok, (1 /)+ Dx, x ir hird2r32) 


E tatoo Cr dador 2 ri3rJ3) 


(2) 


E z k ko,k 
B,(i, ri, ri3)= DM ok Es 


(izo) (3) 


Donde 1, 2 o 3 es el subíndice que 
identifica la presa (La Yesca (1), El Cajón 
(2) o Aguamilpa (3)); A, el producto de las 
probabilidades de transición, q, ,. (j)i=1,3 
en la etapa n, dadas las extracciones ¡A 
y los estados iniciales, 1, y finales, j, de cada 
presa; bpar, la suma del beneficio en la etapa, 
n, dados los estados inicial, i,, i,, 1,, y final j,, 
Jo, j, y SUS respectivas extracciones k,, ss eS 
Bluk2ka (i,,i,,i3), el beneficio esperado hasta 
la etapa n en cada uno de los estados 1, 1,, 1,, 
dadas las extracciones k,, k,, k,; Bara liriaria ), 
el valor esperado óptimo correspondiente a la 
extracción óptima k* en la etapa n + 1. 

Para minimizar el número de cálculos, la 
ecuación (2) se reorganizó como: 
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kay ¡ko yk E 
Br cl = Prot, ko, ke 1112113) 

NS, NS, NS, 
+ y PS > Gx rr h Max, 211291, (13713) 

J¡=1j,=1j3=1 * o 

BaGrrjor3) (4) 
donde: 
NS, 


Pr ko ko ko Ur d2113)= Y dm, ir IA) 
jj=1 
NS» 
+ y mk, (21 )2 mx, rr hi1 142) 


Jal 


NS, 
E DA Qn,k, 37 Ja mxo, rr d11d21 211313 ) (5) 
Ja=1 


Para resolver las ecuaciones (3), (4) y (5), el 
algoritmo de optimización se dividió en dos 
partes, para evitar realizar cálculos repetitivos; 
en la primera sólo se calcula el valor esperado 
del beneficio para cada etapa usando la 
ecuación (5), que se repite de un año a otro, y 
que considera a la función objetivo. Luego, en 
la segunda parte, se supone un valor grande 
N (que puede corresponder a la vida útil del 
sistema); el proceso de cálculo inicia hacia atrás, 
esto es, con la última etapa del año, y en ella 
se suponen los beneficios máximos B* iguales 
a cero. Se aplica la ecuación (4) hasta que la 
diferencia entre la suma de los incrementos 
de los beneficios entre dos años consecutivos 
cumpla con una tolerancia dada que, para fines 
de garantizar convergencia de la solución, se 
supondrá pequeña (10? en este caso). Una vez 
que el proceso converge, se guarda para cada 
presa y para cada etapa la extracción Óptima k* 
y el beneficio total máximo B*. 


Datos de entrada 
Se propuso inicialmente un AV = 200 millones 


de m?, con lo que la capacidad útil de cada presa 
quedó dividida en ocho estados para La Yesca, 


siete para El Cajón y 13 para Aguamilpa. Y en 
un intento por hacer más fina la discretización 
del volumen útil, se consideró un AV = 150 
millones de m', lo que resultó en 10, 9 y 18 
estados para La Yesca, El Cajón y Aguamilpa, 
respectivamente. 

El año se consideró dividido en ocho etapas 
agrupadas en los siguientes meses: etapa uno: 
noviembre; etapa dos: octubre; etapa tres: 
septiembre; etapa cuatro: agosto; etapa cinco: 
julio; etapa seis: junio; etapa siete: febrero- 
marzo-abril-mayo; etapa ocho: diciembre- 
enero. 

Se determinó la probabilidad de los 
volúmenes ingreso para cada una de las 
etapas, a partir de los valores históricos 
registrados en el periodo común de 1981 
a 2008. Para las curvas guía de El Cajón y 
Aguamilpa se utilizaron las elevaciones 
máximas históricas quincenales. Para La 
Yesca se propuso una curva guía similar a la 
de El Cajón, aumentándole un estado, pero 
sin sobrepasar su nivel de aguas máximo 
extraordinario (NAME). Las figuras 3, 4 y 5 
muestran la curva guía para cada presa. 

La curva elevaciones capacidades de cada 
presa, considerando volúmenes a cada AV y 
la elevación con origen en el nivel de aguas 
mínimo ordinario (NAMINO), se representó 
como una ecuación polinomial de segundo 
orden. Así se tiene que: 


Para la Yesca: 
h = -9E-06V? + 0.0529V + 0.2668 (6) 
Para El Cajón: 
h =-7E-06V? + 0.0434V + 0.2391 (7) 
Para Aguamilpa: 
= -7E-07V? + 0.0135V + 0.2546 (8) 
Durante el proceso de optimización se 


consideraron diversos valores para penalizar 
los eventos no deseados en el sistema, que son 
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presencia de derrames, de déficits y de rebases 
de las curvas guía. Se hicieron varios ensayos, 
para cada uno de ellos se obtenía la política 
y con ella se simulaba el funcionamiento del 
sistema. 

Para las simulaciones se consideraron los 
volúmenes quincenales de escurrimiento 
por cuenca propia de las presas La Yesca, 
El Cajón y Aguamilpa (en los cuadros 4a y 
4b se presentan, a manera de resumen, los 
volúmenes mensuales de ingreso consi- 
derados para cada presa). Además, para 
calcular el volumen de evaporación, se 
utilizó el promedio mensual de los valores 
de la lámina de evaporación neta (cuadro 
5). También se consideran datos físicos de 
las presas, como capacidad útil, capacidad 
muerta, capacidad total, nivel de desfogue, 
elevación al NAMINO, volumen máximo 
turbinable, factor de eficiencia, curvas 
elevaciones-capacidades-áreas, etcétera. 

Para simular el sistema con un registro 
mayor que el histórico se generaron diez series 
sintéticas de cien años cada una, usando el 


método de Svanidze modificado (Svanidze, 


1980; Domínguez et al., 2001b; Domínguez et 
al., 2002). 


Resultados y discusión 


Se hicieron diversas combinaciones de valores 
de coeficientes de penalización para los casos 
de déficit, derrame y rebase de la curva guía 
para los dos valores de discretización del 
volumen útil considerados. Para castigar el 
derrame se consideraron valores de 10 000 
y de 100 000; para el déficit, de 100, 1 000 y 
10 000, y para el rebase de la curva guía, de 
1 000. Se hicieron combinaciones de estos 
valores y con cada política que se obtuvo 
se simuló el funcionamiento en forma 
quincenal del sistema. Al analizar el conjunto 
de resultados obtenidos se seleccionaron 
las políticas que mejor desempeño tuvieron 
(maximizar generación de energía y minimizar 
los derrames, déficits y rebases de la curva 
guía). En el cuadro 6 se muestra un resumen de 
los resultados que se obtuvieron al simularlas 
usando el registro histórico. 

De un total de 21 políticas analizadas, las de 
mejor desempeño resultaron la 4, 12, 20 y 21; 
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Cuadro 4a. Ingresos mensuales por cuenca propia en millones de m', La Yesca. 
8 Pp prop: 


Año Enero Feb. Marzo | Abril Mayo | Junio Julio Ago. Sep. Oct. Nov. Dic Total 
1981 60.27 50.99 | 55.64 56.69 82.23 | 198.46 931.66 | 500.21| 434.70| 120.57| 103.03| 135.60| 2730.04 
1982 67.71 50.33 | 46.44 58.87 64.03 65.57 508.73 | 366.77 | 130.33 75.64 57.91 54.83 | 1547.16 
1983 36.843 | 22.432| 17.415| 14.257| 20.312| 297.98 | 1006.843 | 926.713 | 391.307 | 82.416| 41.489 14.31 | 2872.32 
1984 69.87 43.67 | 47.40 45.82 34.09 | 309.19 | 1195.98 | 1107.52 | 596.75| 138.27 76.36 76.40 | 3741.31 
1985 80.85 69.36 | 99.91 85.86 70.43 | 208.80 554.91 | 760.25| 336.98 | 188.28| 113.96 96.38 | 2 665.98 
1986 76.56 60.94 | 64.36 88.81 | 119.92 | 243.34 906.56 | 380.48 | 501.21| 343.65| 156.68| 119.53 | 3062.02 
1987 166.85 | 171.10| 187.21| 146.10| 126.30| 145.73 399.89 | 620.44 | 634.35| 399.17| 184.61| 188.96 | 3370.72 
1988 163.46 | 137.75| 173.62 | 150.39| 156.61| 157.85 737.43 | 1625.02 | 578.61| 164.24 98.59 89.24 | 4232.81 
1989 163.09 81.27 | 64.74 51.11 40.91 49.74 173.07 | 436.12| 296.26 76.82 65.91 68.19 | 1567.21 
1990 126.44| 11223| 81.14 77.20 SIMS6] ELSD:99 572.75 | 2631.54 | 1025.06| 443.06| 154.70 69.40 | 5530.85 
1991 73.91 43.17 | 62.76 68.07 52.89 | 108.28| 4497.60 | 790.97 | 803.19| 220.75 78.62 81.92 | 6882.13 
1992 | 2298.80 | 598.21| 135.50| 136.87| 126.01| 120.71 204.73 | 430.94 | 319.59| 586.88| 129.39| 473.85| 5561.48 
1993 99.30 79.61 | 92.16 74.26 75.14 | 182.16 785.14 | 266.42 | 479.03| 161.71| 105.98 64.16 | 2 465.05 
1994 64.77 54.01 | 78.74| 139.83| 108.34| 186.48 180.32 | 238.63 | 466.72| 163.01 96.34 89.52 | 1866.71 
1995 61.13 58.79 | 89.95 82.17 67.68 | 209.58 643.24 | 1306.44 | 792.58| 195.70| 120.66 96.02 | 3723.93 
1996 58.73 49.08 | 52.07 40.23| 11941| 144.25 257.27 | 268.14 | 618.41 | 537.59 S209 72.13 | 2300.32 
1997 73.35 65.89 | 96.25| 121.38 98.50 | 155.64 314.57| 211.45| 226.72 93.42 82.17 82.09 | 1621.43 
1998 83.39 83.05 | 68.59 33.21 48.37 | 127.59 241.47 | 467.04| 822.25| 620.71| 122.42| 109.87 | 2827.95 
1999 74.388 | 59.419| 62.534| 97.226| 99.239| 158.133 | 444.058 | 375.939 | 439.565 | 143.720| 80.960| 78.280 | 2 113.461 
2000 81.42 78.18 | 67.50 32.60 49.13 | 244.84 171.76| 230.53| 196.63| 128.18 65.09 58.42 | 1404.27 
2001 64.49 20.65 | 55.23 52.38 50.02 | 118.07 561.01 | 436.32 | 550.67 72.49 54.23 46.50 | 2 082.06 
2002 62.12 53.43 | 43.14 34.81 52.34 93.73 527.91 | 452.51| 406.36| 243.05| 138.57 68.85 | 2176.82 
2003 70.15 57.61 | 51.83 39.60 48.11 | 204.98 668.12 | 1262.24 | 1773.77 | 396.76 98.38 60.52 | 4732.06 
2004 90.34 47.89 | 26.35 14.83 18.40 | 451.62 385.29 | 481.77 | 1726.51| 589.85| 102.91 85.25 | 4021.01 
2005 75.59 70.97 | 56.43 38.25 29.05 54.12 279.37 | 438.27 | 379.99| 140.29| 173.23 71.17 | 1806.73 
2006 58.95 591991 IS712 50.38 54071 UA 2297 6901 (83/:821 1887:99; 91.01 793 E2ASiroO 
2007 50.82 28.79 | 38.36 29.37 45.47 | 225.07 931.26 | 1129.77 | 44423| 118.92 65.01 50.93 | 3158.00 
2008 85.72 24.16 | 144.33 44.82 32.79 | 188.66 | 1327.83 | 1723.10| 2292.09| 219.87 78.23 59.62 | 6221.22 
Media 162.12 83.18 | 75.60 68.05 70.42 | 176.01 701.35 | 733.02| 64292| 250.11| 100.69 94.28 | 3157.75 
Son 420.05| 106.11 | 40.89 39.85 36.40 84.67 810.16 | 562.94 | 503.75| 172.89 36.36 80.86 | 1489.51 
pro 5.24 4.55 1.39 0.84 0.66 1.26 4.05 1.81 2.11 1.00 0.69 4.12 1.04 
ri+1,j | 0.9685| 0.4819| 0.7605| 0.8483 | -0.2297| 0.0123 0.2136 | 0.5083 | 0.4025| 0.4699| 0.3987 | -0.0109 
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Cuadro 4b. Ingresos mensuales por cuenca propia en millones de m”, Aguamilpa. 


Año Enero Feb. Marzo | Abril Mayo | Junio Julio Ago. Sep. Oct. Nov. Dic Total 
1981 42.31 16.87 11.69 9.89 5.72 | 105.03| 985.23| 560.41 | 658.23| 135.57 25.40 24.93 | 2581.28 
1982 15.72 10.30 8.19 5.60 4.95 32.17| 468.51 | 318.43 97.92 | 124.27| 103.07| 181.16| 1370.29 
1983 249.65 40.54 46.09 10.23 32.69 57.35| 511.83 | 1059.15| 888.76 | 155.90 61.58 24.80 | 3138.57 
1984 44.98 65.37 12:18 8.14 1152| 218.42 | 1380.50 | 1032.11| 322.20| 101.60 35.96 25.95 | 3258.87 
1985 362.42 34.93 16.70 9.75 6.78| 246.36| 333.05| 843.31| 302.23| 124.02 30.61 50.17 | 2360.34 
1986 6.91 10.66 1.09 0.45 142| 13447| 428.90| 437.39| 504.96 | 239.66 0.10 0.16 | 1766.17 
1987 720.69 | 172.86 85.26 19.41 14.67 47.05| 51297| 841.12| 485.91| 202.86 6.62 39.38 | 3148.81 
1988 21.86 12.87 10.48 8.61 5.65 | 155.71 | 575.09 | 1356.73 | 342.80 Oia 16.25 17.77 | 2602.95 
1989 14.14 9.93 7.39 5.89 4.60 9.27| 237.92 | 612.67 | 371.14 62.84 54.72 51.23 | 1441.74 
1990 17.13 60.75 11.46 5) 8.39| 104.36| 799.58 | 2798.45 | 1039.36 | 365.03 51.49 17.23 | 5278.83 
1991 14.01 10.60 7.46 5.29 3.80 44.32 | 2016.66 | 739.05| 1201.08| 270.12 79.26 98.49 | 4490.14 
119929 [6213391831 [1385:211 86.64 21.24 0.65 30.08 | 251.98| 413.70| 300.61 0.67 57.01 0.11 | 3887.73 
1993 0.50 1.60 0.90 0.98 1.56 0.19| 407.56| 213.52| 804.65 84.88 54.27 6.57 | 1577.18 
1994 28.43 11.89 30.46 0.65 6.26 60.42 | 107.38| 389.37| 340.68| 352.68 27.60 11.08 | 1366.91 
1995 36.15 18.42 40.51 7.54 19.68 78.19 | 2819.99 | 1694.36 | 1057.45 | 147.20 53.76 9.49 | 5982.74 
1996 36.43 39.09 12.69 27.56 43.89 | 178.80| 418.18| 735.35 | 1014.03 | 1135.40 71.65 25.57 | 3738.65 
1997 109.89 65.12 53.79 61.91 68.54 | 221.92| 550.52| 503.70| 497.97| 139.51| 117.72 26.24 | 2416.82 
1998 28.16 23.42 2010217) 18.66 58.36 72.29| 416.17| 824.96| 994.02| 502.88 48.91 21.16 | 3030.24 
1999 11.03 13.11 4.55 8.92 32.16| 175.19| 90244| 624.41| 1090.22| 130.35 10.84 12.19 | 3015.41 
2000 16.81 5:79) 5.46 la 5.63 | 220.30| 233.41| 279.24| 153.35 79.29 19.68 16.42 | 1040.45 
2001 15.90 4.88 4.82 5.18 5.40 54.80 | 376.00| 463.71| 420.20 8.36 15.83 16.83 | 1391.91 
2002 8.81 32.78 27.82 0.13 7.37 | 102.21| 718.76 | 1085.43 | 670.89 | 405.38 79.56 18.75 | 3157.88 
2003 78.23 7.03 6.61 7.75 3.39| 115.49 | 636.28 | 1636.68 | 397.72| 435.09 57.41 14.17 | 3395.85 
2004 79.69 9.82 0.38 8.95 57.50 | 400.87 | 570.37 | 1180.83 | 2698.27 | 338.06 5/39 37.99 | 5440.09 
2005 43.48 | 113.03 0.30 0.65 17.39 20.26| 503.08| 918.47 | 516.40| 139.32 36.78 26.24 | 2335.39 
2006 8.24 43.00 7.01 16.83 8.43 SiS] AS i914 | ESSO51S 252.031 ISS 36.78 26.24 | 1194.18 
2007 22.51 12.90 12.69 4.38| 17997| 11455| 534.67 | 572.14| 468.49 85.87 36.78 26.24 | 2071.20 
2008 78.23 16.31 12.69 8.95 15.37| 132.26| 860.25 | 1395.83 | 613.26| 126.61 22.06 55.94 | 3337.75 
Media 159.01 44.61 19.52 10.51 22.56| 11455| 668.08| 853.20| 660.89 | 220.89 45.32 31.52 | 2850.7 
Desv. 
cd 452.29 76.39 23.22 12.10 36.17 89.09 | 574.62 | 555.11| 505.79| 221.66 28.19 35.41 | 1315.9 
pe 4.53 3.70 1.94 3.09 3.37 1.34 2.55 1.78 2.54 2.81 0.70 3.19 0.74 
rj+1, 0.9375 | 0.7443| 0.4912 | 0.2574| 0.2472 | -0.0009| 0.4089| 0.3256| 0.3369| 0.2522| 0.4386| 0.3896 
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Cuadro 5. Evaporación mensual neta (mm). 


Presa Ene. Feb. Mar. Abr. May. Jun. Jul. Agos. | Sept. Oct. Nov. Dic. 
La Yesca 69.59 | 104.71 | 179.03 | 211.42 | 223.49 | 31.27 |-133.32 | -96.67 | -49.58 | 48.85 | 68.79 | 50.804 
El Cajón 69 95 143 170 191 5% 11215 -103 500) 43 83 54 

Aguamilpa | 23.5 63.1 97.4 115.6 115 -40.6 | -269.1 | -221.1 | -154.6 0.1 33.2 24.3 
Cuadro 6. Resumen de la simulación conjunta del registro histórico. 
Sistema hidroeléctrico del río Santiago. Periodo: 1981-2008. 
Política | Energía promedio anual GWh E Derrame, 105 m' Déficit, 10% de m' 
La Yesca | El Cajón | Aguamilpa | GWh | La Yesca | El Cajón | Aguamilpa | Total | La Yesca | El Cajón | Aguamilpa | Total 
4 1 145.5 1 101.6 2 020.8 4267.9 | 5442.0 | 2222.1 4 297.7 119618 | 300.5 0.0 62.3 362.8 
9 1 064.2 1116.7 2 002.8 4183.7 | 3871.0 | 2035.3 2.295.9 8202.1 | 17066.9 | 1391.4 788.4 19 246.7 
12 1 136.4 1 143.1 2 045.0 4324.6 | 5170.6 | 2464.9 5:528.7 13 164.3 38.2 0.0 0.0 38.2 
Limitando la curva guía 
20 1 129.2 1 084.1 2 007.1 4220.4 | 50322 | 2309.2 4 472.2 11813.5 | 965.7 374.6 140.1 1 480.3 
21 1 136.4 1 143.1 2 045.0 4324.6 | 5170.6 | 2464.9 5: 528.7 13 164.3 38.2 0.0 0.0 38.2 
Cuadro 7. Coeficientes de penalización considerados para cada presa. 
Coef. derrame Coef. déficit Coef. rebase curva guía 
Presa Política Política Política 
4 9 4 9 4 
1 10 000 100 000 10 000 1000 1000 1000 
2 10 000 100 000 100 100 1000 1000 
3 10 000 100 000 100 100 1000 1000 


en todas ellas, los coeficientes de penalización 
fueron los mostrados en el cuadro 7 para la 
política 4. 

Se incluyó también en el cuadro 6 la política 
9, que resulta al imponer las penalizacio- 
nes indicadas en el cuadro 7, para mostrar 
cómo cambian los resultados al modificar las 
penalizaciones. Así, al comparar los resultados 
de las políticas 4 y 9, se observa que al aumentar 
la penalización por derrame y disminuir la de 
déficit en la presa 1, y aumentar la de derrame 
en las presas 2 y 3, la energía total generada 
disminuye muy ligeramente, mientras que el 
derrame total en el sistema disminuye un 30%, 
pero el déficit se incrementa exageradamente, a 


niveles que la CFE consideró inaceptables, por 
lo que la política 9 se descartó para los análisis 
posteriores. 

Las políticas 4 y 9 se determinaron usando 
una discretización de la capacidad útil de las 
presas de 200 millones de m* y la 12 con un 
AV de 150 millones de mY. Las políticas 20 y 21 
corresponden a las discretizaciones de la 4 y la 
12, respectivamente, pero en ambas se limitó la 
curva guía de la presa Aguamilpa de manera 
que no sobrepasara su NAMO en ninguna 
quincena. 

De acuerdo con los resultados del cuadro 6, 
la reducción en el tamaño de la discretización 
no conduce a cambios cualitativos en los 
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resultados. Con la política 4 se presenta un 
derrame total cercano a los 12 000 millones 
de m”, mientras que con la 12 resultan poco 
más de 13 000 millones de m*; en cambio, con 
la política 4, se presenta un déficit mayor que 
con la 12. De aquí se puede concluir que si lo 
que se pretende es controlar déficits pagando 
con presencia de derrames, la política 12 es 
mejor que la 4; sin embargo, si es lo contrario, 
la política 4 derrama mil millones menos que 
la 12, aunque no puede evitar que las presas 
de La Yesca y Aguamilpa lleguen a su nivel 
mínimo y presenten déficit. 

Al limitar la curva guía de la presa Agua- 
milpa con el NAMO, se observan cambios 
en la política 20 (respecto a la 4); como era 
de esperarse, la política 20 controla mejor 
los derrames, pero paga el precio con un 
incremento en el déficit en las tres presas. La 
política 21 se mantiene sin cambio respecto a 
la 12. Se podría decir con esto que el limitar la 
curva guía de la presa Aguamilpa no represen- 
ta un cambio significativo en los resultados; 
sin embargo, el hecho de que se limite la curva 
guía al nivel del NAMO es más congruente con 
las condiciones del diseño original del embalse 
y hace más seguro el funcionamiento de la obra 
de excedencias. 

Respecto a la generación de energía pro- 
medio anual, con objeto de tener un parámetro 
cualitativo de comparación, se utilizó la 
estimación reportada por la CFE (CFE, 2010), 
que para El Cajón es de 1 228 GWh, mientras 
que las simulaciones con las políticas 4 y 12 
reportan 1 102 y 1 143 GWh, respectivamente; 
es decir, que se tiene una diferencia de 126 GWh 
con la política 4, la que se reduce a 85 GWh con 
la política 12. Para la presa de Aguamilpa, la 
generación media anual estimada por la CFE 
es de 2 133 GWh y los resultados que alcanzan 
las dos mejores políticas (la 4 y la 12) es de 2 021 
y 2 045 GWh, lo que da una diferencia de 112 y 
88 GWh, respectivamente. La comparación con 
las estimaciones de la CFE es sólo cualitativa 
debido a que en la referencia utilizada no 
se indican los datos de escurrimiento y 


evaporación ni la eficiencia de generación que 
utilizaron. 

Usando el registro histórico (28 años) para 
simular el sistema, en la figura 6 se muestra 
que con la política 20, en La Yesca casi cada 
cinco años ocurren derrames, cada nueve años 
en El Cajón y casi cada seis años en Aguamilpa. 
La política 21 disminuye la frecuencia de 
derrames en la primera presa a casi cada seis 
años; la mantiene igual en El Cajón (casi cada 
nueve años), y en Aguamilpa la incrementa un 
poco, casi cada cinco años. De la misma figura 
6 se tiene que la ocurrencia de déficit para la 
política 20 es de cada nueve años para La Yesca, 
no se presentan en El Cajón y ocurren cada 28 
años en Aguamilpa. La política 21 logra que 
en El Cajón y Aguamilpa no se tengan déficits, 
y que cada 28 años se presente el evento en 
La Yesca. Los rebases de la curva guía para la 
política 20 se presentan en promedio casi cada 
cinco años en La Yesca y cada 3.5 y 4 años 
en El Cajón y Aguamilpa, respectivamente. 
Con la política 21, La Yesca se mantiene sin 
cambio, pero El Cajón y Aguamilpa presentan 
ahora prácticamente cada dos años rebase de 
los límites de la curva guía. 

Para tener un periodo más confiable de 
simulación del sistema, que aporte mayor 
detalle respecto a los posibles eventos extre- 
mos (situaciones de derrames y déficits), se 
generaron diez series sintéticas de cien años 
cada una, y se hizo con ellas la simulación 
de las políticas 20 y 21. Al comparar los 
resultados de las simulaciones con la política 
20 (cuadro 8) con los de la política 21 (cuadro 
9), se tendría con la primera un derrame 
promedio de 35 213.5 millones de m?* en cien 
años para la suma de las tres presas, contra 39 
865.2 millones de mY, para la política 21; esto 
es, un 13% de aumento en la magnitud del 
derrame. En el caso del déficit, el promedio 
en el sistema sería de 10 258.3 millones de 
m' en cien años con la política 20, contra los 
1 100.1 millones de 1m* de la 21; en este caso, 
el déficit disminuiría en cerca del 90%; es 
decir, con la política 21 se tendría un déficit 
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Número de años con derrame, déficit y rebase de CG (registro histórico) 
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Figura 6. Número de años con derrames, déficits y rebases de la curva guía para las políticas 20 y 21, 
usando el registro histórico. 


Cuadro 8. Resumen de la simulación del funcionamiento del sistema con la política 20, y diez series sintéticas de 


cien años cada una, río Santiago, Nayarit. 


Política 20 Energía promedio anual (GWh) Derrame total, 10* m* Déficit total, 10% de m* 
Serie La Yesca |El Cajón | Aguamilpa | Total | La Yesca | El Cajón | Aguamilpa | Total | La Yesca | El Cajón | Aguamilpa | Total 

dl 1 092.2 1016.2 1 857.6 3 966.0 3 696.4 431.3 207936 | 24921.3 | 5137.7 2ATOLO 1709.3 9 598.9 
2 1139.5 1 057.8 19771 41744 | 13858.9 | 6953.0 22 613.7 | 43425.6 | 3838.4 979.6 977.2 5795.2 
3 1 145.8 1 088.2 19195 41474 | 148045 | 8133.6 253295 [|[48267.5 | 8 109.9 3 895.0 1821.6 13 826.5 
4 1140.2 1074.2 1 918.8 4 133.3 93472 | 33424 111821 |23871.7 | 60946 2 973.9 1 404.9 10 473.4 
E 1112.6 1 052.4 1923.6 4 088.6 7 461.1 1005.5 10650.0 | 19116.6 | 3329.1 934.9 906.8 5 170.8 
6 978.5 931.9 1799.32 3709.7 | 12828.9 | 4136.0 19001.5 | 35966.4 | 17 835.4 8953.7 | 2005.1 28 794.2 
dl 1 209.4 1 134.2 2017.7 4361.3 | 13943.6 | 9036.6 34 908.7 | 57888.9 | 4830.2 1 457.0 1 032.4 7319.7 
8 1162.1 1 108.3 1960.1 4230.5 | 13395.0 | 4746.9 31429.6 | 495715 | 4439.1 16144 1263.6 7317.1 
Ls 1 080.2 1011.8 1 853.8 3945.8 | 114411 | 72355 17 095.0 | 35771.6 | 4685.9 1 636.6 51.5 6.374.0 
10 1 092.2 1015.9 1976.4 4 084.6 3 696.4 431.3 9206.7 | 133344 | 51523 2 329.1 432.5 7913.9 

Promedio | 1115.3 1049.1 1 919.8 40842 | 104473 | 4545.2 202210 [352135 | 63452 2 752.6 1160.5 10 258.3 
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Cuadro 9. Resumen de la simulación del funcionamiento del sistema con la política 21, y diez series sintéticas de 


cien años cada una, río Santiago, Nayarit. 


Política 21 Energía promedio anual (GWh) Derrame total, 10% m* Déficit total, 10* de m* 
Serie La Yesca | El Cajón | Aguamilpa | Total | La Yesca | El Cajón | Aguamilpa | Total La Yesca | El Cajón | Aguamilpa | Total 

1 1102.6 | 1076.9 1909.7 4 083.1 39914 | 1023.6 | 22248.6 | 27 263.7 289.7 0.0 0.0 289.7 
2 1148.9 | 1129.0 2 017.0 42948 | 140123 | 7750.9 | 25632.1 47 395.3 292.5 0.0 0.0 292.5 
3 11153:9 | 11366 11956.2 4246.8 | 15223.2 | 10418.1 282690 | 539103 | 19761 0.0 0.0 1976.1 
4 11448 | 1129.2 1 965.6 4 239.6 9757.8 | 43722 13 175.1 27305.1 | 1280.1 0.0 4.7 1 284.8 
iS 11194 | 1104.5 1961.8 4 185.6 8054.1 | 3089.0 15551.3 | 26694.4 683.8 0.0 s2.7 716.5 
6 987.8 988.8 1836.5 3813.1 | 13230.0 | 62264 | 237224 | 43178.8 | 4559.6 0.0 5.6 4 565.2 
J 1213.0 | 1188.5 2. 064.7 4466.2 | 14945.4 | 10310.2 | 36965.1 62 220.7 551.2 0.0 0.0 S51.2 
8 1170.7 | 11575 1 995.4 4323.6 | 13709.8 | 6687.9 | 35414.8 55 812.5 466.6 0.0 0.0 466.6 
e 1088.4 | 1073.3 1 902.2 4063.9 | 11436.8 | 7613.1 18 942.6 37 992.5 560.2 0.0 0.0 560.2 
10 1102.1 | 10735 20225 4 198.1 4193.8 | 11675 11 517.8 16 879.1 297.8 0.0 0.0 297.8 

Promedio | 1123.2 | 1105.8 1962.6 4191.5 | 10855.4 | 5865.9 | 23143.9 39865.2 | 1095.8 0.0 4.3 1100.1 


mucho menor en el sistema. Es interesante 
resaltar que estos resultados son congruentes 
con los obtenidos para el registro histórico, 
lo que confirma la validez de las secuencias 
generadas. 

Las curvas guía diseñadas se rebasan 
en varias ocasiones durante la simulación 
del funcionamiento del sistema con ambas 
políticas. La figura 7 muestra el promedio 
en los cien años de las diez series sintéticas 
generadas. 

Se puede apreciar que la política 20 tiene un 
mejor control sobre este evento que la política 
21; con la primera se rebasa la curva guía en La 
Yesca cada 1.7 años, 2.6 años en El Cajón y casi 
cada, 1.4 años, en Aguamilpa; mientras que 
con la política 21, en La Yesca, cada 1.6 años se 
rebasa la curva guía; cada 3.4 años en El Cajón, 
y casi cada cinco años en Aguamilpa. 

En cuanto a generación de energía, con la 
política 20, las dos presas (El Cajón y Aguamil- 
pa) generan menos que el valor estimado por la 
CFE; a la primera le falta aproximadamente un 
17% y a la segunda un 11%. La política 21 tiene 
un mejor desempeño respecto a la generación 
de energía, pero tampoco alcanza los valores 
históricos estimados por la CFE, se queda con 


un faltante de 11% para El Cajón y un 8.6% 
para Aguamilpa. 

De acuerdo con los cuadros 6, 8 y 9, la 
discretización en más estados, esto es, pasar 
de un AV de 200 a uno de 150 millones de 
n*Y, conduce a resultados más precisos, pero 
cualitativamente similares; sin embargo, en 
los recursos computacionales usados existen 
diferencias marcadas, principalmente en 
los tamaños de los arreglos —mismos que 
condicionan al compilador al momento de 
crear el programa ejecutable— y en el tiempo 
de ejecución de los programas de cómputo 
empleados, que pasó de dos minutos (120 s) 
para el AV de 200 millones de m' a 10 minutos 
(600 s) para la discretización más fina de 150 
millones de mí; es decir, que al disminuir en un 
25% el tamaño de la discretización, el tiempo 
de ejecución se multiplicó por cinco. Debido a 
ello, no se consideró necesario usar intervalos 


de discretización menores. 
Conclusiones 
Se logró desarrollar un algoritmo para re- 


solver el sistema de ecuaciones resultante de 
la aplicación de la programación dinámica 
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Promedio de años con rebase de las curvas guía (series sintéticas) 
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Figura 7. Promedio de número de años en los que se rebasa las curvas guía en el sistema para las diez series sintéticas. 


estocástica para un sistema de tres presas en 
cascada que funcionara con éxito, a pesar 
del gran número de estados que implica (el 
número de estados necesarios para carac- 
terizar el sistema es igual al producto de los 
estados correspondiente a cada presa). Con 
el algoritmo desarrollado se determinaron 
políticas de operación de un sistema de presas 
ubicado en el río Santiago, con la propuesta de 
dos incrementos de volumen (AV = 200 y AV= 
150 millones de m?*) para dividir la capacidad 
útil de las presas y considerando curvas guía 
obtenidas a partir de la operación histórica. 

Las políticas de operación obtenidas con- 
cilian el objetivo de la generación de energía 
con el de evitar derrames o déficit, tomando en 
cuenta el efecto de las curvas guía. 

En este estudio, el tiempo de cómputo 
empleado no era un factor muy importante a 
considerar; se le dio más peso a que el algoritmo 
propuesto usara en forma óptima la estructura 
misma de la programación dinámica, al dividir 
el problema en subproblemas más pequeños, 
resolverlos usando la recursividad y usar las 
soluciones para obtener, por el mismo principio 
del método, una solución óptima al problema 
completo. 

En el cuadro 10 se presenta el ejemplo de 
la política de operación 4, para la etapa de 


noviembre, en la que se indica que (por razones 
de espacio, en el cuadro sólo se muestran los 
primeros tres estados) si la presa uno está en 
el estado 01, la presa dos está en el estado 5 
(renglón 5) y la presa tres está en el estado 6 
(columna 6); la política indica extraer 010201, 
que significa extraer en la quincena un volumen 
01 x 100 = 100 millones de m* de la presa La 
Yesca, 02 x 100 = 200 millones de m' de la presa 
El Cajón y 01 x 100 = 100 millones de m', de la 
presa Aguamilpa. 
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Cuadro 10. Política de operación 4 para la etapa 1, noviembre. Cada unidad corresponde a una extracción 


de 100 millones de m* por quincena. 


Etapa 1: noviembre 


Presa: 1, estado: 01 


010101 


010102 
010101 
010101 
010101 
010201 
010201 
010301 


010103 
010101 
010101 
010101 
010201 
010201 
010401 


010104 
010101 
010101 
010101 
010201 
010201 
010401 


010105 
010101 
010101 
010101 
010201 
010201 
010401 


010106 
010101 
010101 
010101 
010201 
010301 
010401 


010107 
010101 
010101 
010101 
010201 
010301 
010401 


010108 
010101 
010101 
010101 
010201 
010301 
010401 


010109 
010101 
010101 
010101 
010201 
010301 
010401 


010101 
010101 
010101 
010101 
010201 
010301 
010402 


010101 
010101 
010101 
010101 
010201 
010302 
010403 


010101 
010101 
010101 
010101 
010202 
010303 
010404 


010102 
010102 
010102 
010102 
010203 
010304 
010405 


010101 
010101 
010101 
010201 
010201 
010301 


XN[DO|OI|2|O0|Nn|Re 


Presa: 1, estado: 02 


1 2 3 4 5 6 7 8 9 10 11 12 13 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 [ 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010201 | 010201 | 010201 | 010101 | 010203 
010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010302 | 010303 | 010304 
010201 | 010201 | 010201 | 010201 | 010201 | 010401 | 010401 | 010401 | 010401 | 010402 | 010403 | 010404 | 010405 
010301 | 010301 | 010501 | 010501 | 010501 | 010501 | 010501 | 010501 | 010502 | 010503 | 010504 | 010505 | 010506 


XN|[DO|OIO|[2R|O0OQ|nNn]|Re 


Presa: 1, estado: 03 


1 2 3 4 S 6 Y 8 9 10 11 12 13 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 ¡ 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010101 | 010102 
010201 | 010201 | 010201 | 010201 | 010301 | 010301 | 010301 | 010301 | 010301 | 010301 | 010302 | 010303 | 010304 
010201 | 010201 | 010201 | 010401 | 010401 | 010401 | 010401 | 010501 | 010401 | 010402 | 010403 | 010404 | 010405 
010501 | 010501 | 010501 | 010501 | 010601 | 010601 | 010601 | 010501 | 010502 | 010503 | 010504 | 010606 | 010607 


XN[DO|O|(2|O0|nNn]| ea 
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Abstract 


ARGANIS-JUÁREZ, M.L., MENDOZA-RAMÍREZ, R. € DOMÍNGUEZ-MORA, 
R. Operation of three hydroelectric dams using guide curves and stochastic dynamic 
programming. Water Technology and Sciences (in Spanish). Vol. IL, No. 3, July- 
September, 2012, pp. 97-114. 


Dynamic Programming (DP) is an optimization procedure developed in the 1950s by 
Richard Bellman (Bellman, 1957). It is a mathematical tool that has been applied in several 
fields, including engineering, artificial intelligence and economics, among others (Dormido 
et al., 2002). The area of water resources planning is no exception and applications with 
good results have been reported in several countries. A limitation is high computing costs, 
the need for many resources both in terms of time as well as memory, with the demand of 
these resources growing exponentially as the number of state variables increases. This work 
presents a way in which to organize a computing program to handle this problem as well as 
to determine policies for the optimal operations of a system of three serial hydropower dams 
by varying the penalties for deficits, spills and exceeding guide curves (the storage levels not 
to be surpassed during operation for the safety of the system, which are established by the 
National Water Commission). The policies obtained were simulated using historical records 
and synthetic records to test their benefit to the operations of the system. 


Keywords: dynamic programming, serial operation, operations policies, guide curves. 
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